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1. Introduction 

1.1. Historical perspective 

Strange objects, which are poorly characterised by their topological dimension, such as 
the Weierstrass continuous, but nowhere different iable function, the van Koch curve, 
the Sierpihski gasket (see Figured]), the Sierpihski carpet (see Figure H]), etc., have been 
familiar in mathematics for a long time. Now such objects are called fractals. 

The word "fractal" was coined by Mandelbrot in the 1970s. His foundational 
treatise [1] contains a great diversity of examples from mathematics and natural sciences. 
Coastlines, topographical surfaces, turbulence in fluid, and so on are just few instances 
from the variety of natural objects that may be described as fractals. There is no 
generally agreed exact definition of the word "fractal" . 

The fractals studied in the context of analysis on fractals are all self-similar. 
Moreover, we deal here mainly with finitely ramified fractals, i.e., fractals that can 
be disconnected by removing a specific finite number of points. 

Initial interest in processes (analysis) on fractals came from physicists working in 
the theory of disordered media. It turns out that heat and wave transfer in disordered 
media (such as polymers, fractured and porous rocks, amorphous semiconductors, etc.) 
can be adequately modelled by means of fractals and random walks on them. (See the 
initial papers by Alexander and Orbach [2], Rammal and Toulouse [3]. See also the 
survey by Havlin and Ben-Avraham [1| and the book by the same authors [5] - for an 
overview of the now very substantial physics literature and bibliography.) 

Motivated by these works, mathematicians got interested in developing the of 
"analysis on fractals" . For instance, in order to analyse how heat diffuses in a material 
with fractal structure, one needs to define a "heat equation" and a "Laplacian" on a 
fractal. The problem contains somewhat contradictory factors. Indeed, fractals like the 
Sierpihski gasket, or the van Koch curve do not have any smooth structures and one 
cannot define differential operators on them directly. We will comment on the difference 
between the euclidean and the fractal situation throughout this paper. 

1.2. Probabilistic approach 

In the mid 1980s probabilists solved the problem by constructing "Brownian motion" on 
the Sierpihski gasket. Goldstein [6], Kusuoka [7], and a bit later Barlow and Perkins [8], 
independently took the first step in the mathematical development of the problem. Their 
method of construction is now called the probabilistic approach. Namely, they considered 
a sequence of random walks X^"^ on graphs G„, which approximate the Sierpihski gasket 
F, and showed that by taking a certain scaling factor, those random walks converge to 
a diffusion process on the Sierpihski gasket. 

In order to obtain a nontrivial limit, the time scale for each step should be the 
average time for the random walk on Gn+i, starting from a point in to arrive at a 
point in Gn, except for the starting point. By the self-similarity and symmetry of F this 
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average time r is independent of n and it is equal to the average time for the random 
walk on Gi, starting from a to arrive at either b or c. (Here a, b, c, are the vertices of 
the initial equilateral triangle Go-) Elementary calculations show that r = 5. According 
to [SllTllH] the processes 2~^X^'^\[5^t]) weakly converge (as n — )• oo) to a non-trivial 
limit Xt on F, which is called Brownian motion on F. In this approach the Laplacian is 
defined as an infinitesimal generator of Xt. 

An important observation was made by Barlow and Perkins in [8]. Let Z„ = T" *^ 
be the first hit time by X^""^ on Gq. Then Zn is a simple branching process. Its off-spring 
distribution rj has the generating function q{z) = E(2'') = z'^/ {A — 3z) and, in particular, 
E(77) = ^'(1) = 5. Thus, Zn is a super- critical branching process. It is known (see [9]) 
that in this case 5~"Z„ tends to a limiting random variable Z^q. 

The moment generating function of this random variable 

f{z) = Ee-^^°° 

satisfies the functional equation 

f{Xz) = qUXz)), (1) 

which is the Poincare equation (see also Section [5] below). 

The Poincare equation associated with the Brownian motion turns out to be a 
very useful tool in the study of the detailed properties (for example, heat kernel) of the 
Brownian motion. 

Lindstr0m [10] extended the construction of the Brownian motion from the 
Sierpihski gasket to more general nested fractals., (which are finitely ramified self-similar 
fractals with strong symmetry). The Lindstr0m snowflake is a typical example of a 
nested fractal (see Figure [2]). 

Readers may refer to Barlow's lecture notes [llj for a self-contained survey of the 
probabilistic approach. 

1.3. Anomalous diffusion 

It has been discovered in an early stage already (see [21 El El O [8]) that diffusion on 
fractals is anomalous, different than that in a regular space. For a regular diffusion, 
or (equivalently) a simple random walk in all integer dimensions d, mean-square 
displacement is proportional to the number of steps n: E(X„)^ = en (Pick's law, 1855). 
On the other hand, in the case of the Sierpihski gasket W{Xn — x)^ x ra^/^, where x 
means that the ratio between the two sides is 'bounded above and below by positive 
constants' and /3 = Ig5/lg2 is called the walk dimension. 

This slowing down of the diffusion is caused, roughly speaking, by the removal of 
large parts of the space. 

De Gennes [12] was amongst the first, who realised the broad importance of 
anomalous diffusion and coined the suggestive term "the ant in the labyrinth" , describing 
the meandering of a random walker in percolation clusters. 
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1.4- Analytic approach 

The second approach, based on difference operators, is due to Kigami [13]. Instead of 
the sequence of random walks, one can consider a sequence of discrete Laplacians on 
a sequence of graphs, approximating the fractal. It is possible to prove that under a 
proper scaling these discrete Laplacians would converge to an "well-behaved" operator 
with dense domain, called the Laplacian on the Sierpihski gasket. This alternative 
approach is usually called the analytic approach. 

Later it was extended by Kigami p ^ [T5 | [T6] to more general class of fractals - post 
critically finite self-similar sets (p.c.f), which roughly correspond to finitely ramified 
self-similar fractals. 

The two approaches described above are complementary to each other. 

The advantage of the analytic approach is that one gets concrete and direct 
description of harmonic functions, Laplacians, Dirichlet forms, etc. (See also pTtlTS].) 

On the other hand, however, the probabilistic approach is better suited for the 
study of heat kernels. Moreover, this approach can be applied to infinitely ramified 
self-similar fractals, which include the Sierpihski carpet, as a typical example (cf., [T9]). 

1.5. Functional equations in the analysis on fractals 

One important example, where the Poincare equation arises in connection with analysis 
on fractals, has been mentioned already at the end of Section 11.21 

Further applications of functional equations in this field are related to spectral zeta 
function of the Sierpihski gasket and other more general fractals having spectral 
decimation. The phenomenon of spectral decimation was first observed and studied by 
Fukushima and Shima ( [T7ll20l[2T]) and further progress has been made by Malozemov 
and Teplyaev [22] and Strichartz [23]. 

The definition of spectral decimation is given in Section 14.21 (see Definition [8] 
below). It implies, in particular, that eigenvalues of the Laplacian A on the fractal, 
which admits spectral decimation, can be calculated by means of a certain polynomial 
p{z), or rational function R{z). Hence, the spectral zeta function may be defined by 
means of iterations p'-"^ of p, or i?*^"^ of R (see [21] )• 

The above mentioned iteration process, as is well known in iteration theory, may 
be conveniently described by the corresponding Poincare equation: 



where X = p'{0) > 1 (see, for example, Beardon [25] or Milnor |26j). 



Using that, in Section 14.41 we obtain the meromorphic continuation of the zeta 
function (/^ into the whole complex plane on the basis of the knowledge of the asymptotic 
behaviour of the Poincare function $ in certain angular regions. 

The poles of the spectral zeta function are called the spectral dimensions (see 
[271 [281 [29]). Foi' the physical consequences of complex dimensions of fractals - see 




(2) 
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( [301 EI])- In Section I4.5.H we use the Poincare function $ for the calculation of 
Casimir energy on a fractal. 

Finally, one can expect that functional equations with rescaling naturally come 
about from problems, where renormalisation type arguments are used to study self- 
similarity. Furthermore, functional equations, in contrast with differential ones, do not 
require any smoothness of solutions; they may possess nowhere differential solutions, for 
example. 

1.6. Notes and remarks 

There is now a number of excellent books, lecture notes and surveys on different aspects 
of analysis and probability on fractals [5l[IIl[l6l[29l|32l|33l[3ai35l[36]. 

It is next to impossible to describe all activities in this area. The objective of the 
present paper is different. We restrict ourselves to a brief overview of various approaches 
in the study of the Laplacian and its spectral properties on certain self-similar fractals, 
and we put much emphasis on the deep connection between the latter problem and 
functional equations with rescaling, and the classical Poincare equation, in particular. 

The Poincare equation plays a very important role in the mathematical theory of 
dynamical systems and in iteration theory, in particular, but is still much less known 
in the physics literature. Hopefully, the current presentation of this realm of problems 
intended for a general audience may fill this gap. 

We begin our presentation of the Poincare equation in Section [5] with a brief 
description of the general case, when the coefficients of the polynomial P and the scaling 
factor A are complex numbers, and only afterwards turn to a detailed discussion of the 
real case. So far, the real case only has arisen in analysis on fractals. However, we think 
that the general case is interesting on its own, and, probably, will also find applications 
in the future. 

In this overview we do not or only cursorily touch the following important topics: 

• analysis on infinitely ramified fractals such as the Sierpiski carpet [19l[37]. The 
very recent progress that has been made in proving the uniqueness of the diffusion 
on the Sierpihski carpets |38] provided a unification of the different approaches to 
diffusion on this class of fractals. 

• heat kernel long time behaviour and Harnack inequalities on general underlying 
spaces, including Riemannian manifolds, graphs and fractals, as special cases 
[39111011111112] 

• analysis by means of potential theory and functional spaces (Sobolev or Besov 
spaces) techniques [33lll3llllll5lll6lll71ll8] 

We refer the reader to the original papers. 
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2. Fractals and iterated function systems 

Let us first recall that for any finite set of linear contractions on M'^ 

i^j(x) = bi + Ai(x - bi), i = l,...,m 

with fixed points bj and contraction matrices Aj (|| Aj|| < 1, i = 1, . . . , m) there exists 
a unique compact set K cM.'^ satisfying 

m 

K = [jF,{K). (3) 

i=l 

In general, the set K obtained as fixed point in ([3]) is a "fractal" set (cf. [l9l[50]). If the 
matrices Aj are only assumed to be contracting, finding the Hausdorff-dimension of K 
is a rather delicate (cf. [SI1E2]) problem. 

For this paper we will make the following additional assumptions on the family of 
contractions F = {Fj \ i = 1, . . . , m}: 

(i) the matrices Aj are similitudes with factors < 1: 

Vx G M'^ : ||Aix|| = ai||x|| 

(ii) F satisfies the open-set-condition (cf. namely, there exists a bounded open 
set O such that 

m 

\Jf,{o)co 

i=l 

with the union disjoint. 

Assuming that the maps in F are similitudes and satisfy the open-set-condition, 
the Hausdorff-dimension of the compact set K given by equals the unique positive 
solution s = p of the equation 



m 



J2^t = l- (4) 



i=l 



The projection map 

vr: {l,...,mf 

(ei, £2, • • •) ^ lim„^oo o o ■ ■ ■ o F,^ (x) 

(the limit is independent of x G M.'^) defines a "parametrisation" of K. If {1, . . . ,m}^ 
is endowed with the infinite product measure fi given by 

/i ({(ei, 62, . . .) \ El = ii,e2 = i2,- ■ - Sk = ik}) = (anttia " " " J'' > (6) 

then TT is a //-almost sure bijection. The normalised p-dimensional Hausdorff-measure 
on K is then given by 

-H'KiE) := ^^^J§^ = -MiE) = P {n^-'\E)) . (7) 
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This is the unique normahsed measure satisfying (cf. [53^^ Theorem 28]) 

m 

'n^{E) = Y.<'^K{Ft"\E)). (8) 

1=1 

By the above discussion the functions 

5„(x) = (7r(-i)(x))^ 

are defined 7{^-almost everywhere on K. We set 

5„(x)=min(7r(-i)({x})),„ (9) 

to define them everywhere. The set K together with the address space S = {1, . . . , m}^ 
and the maps Fi (i = 1, . . . , m) defines a self-similar structure {K, S, (-Fi)^i). 

In order to allow sensible analysis on the fractal set K, we need some further 
properties. Especially, since we will later study diffusion processes, we need K to 
be connected. On the other hand, the techniques introduced later require a finite 
ramification property usually called post-critical finiteness (p. c. f.). 

Connectivity of K is characterised by the fact that for any pair there exist 

ii, 12, . . . ,in £ {1, • • • 5 ^} with i = ii and j = in such that Fi^[K) fl Fi^^^{K) ^ for 
£= (cf. [54J). 

Definition 1 Let {K,J], be a self-similar structure. Then the set 



is called the critical set of K. The post-critical set of K is defined by 

oo 
n=l 

where a : S — t- S denotes the shift map on the address space S. If P is a finite set, 
then [K,!], {Fi)^i) is called post-critically finite (p. c. f.). This is equivalent to the 
finiteness of C together with the fact that all points of C are ultimately periodic. 

The following sequence of finite sets will be used in Section 13.51 to define a 
sequence of electrical networks giving a harmonic structure on K. For more details we 
refer to [T6l, Chapter 1] 

Definition 2 Let {K, S, (-Fi)™ i) be a post- critically finite self-similar structure and P 
its post-critical set. Let Vq = vr(P) and define Vm iteratively by 



Vn+l = [jF,{Vn). 



i=l 



The sets Vn are then finite, increasing (Vn C Vn+i) and 
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Figure 1. The Sierpinski gasket; the points of Vq are the vertices of the triangle 



Figure 2. The Lindstr0m snowflake with the corresponding set Vq 



3. Laplace operators on fractals 

3.1. Laplace operators on compact manifolds 

Before we introduce the Laplace operator on certain classes of self-similar fractals, let 
us shortly discuss the situation in the manifold case, because this gives the motivation 
for the different approaches in the case of fractals. Let M be a compact Riemannian 
manifold with a Riemannian metric g given as a quadratic form Qx on the tangent space 
TxM for X G M. As usual, we assume that the dependence of on x is diffcrentiable. 
Then the quadratic form g^ defines an isomorphism ag between the tangent space T^M 
and its dual T*M (and thus on the tangent bundle TM and the cotangent bundle 
T*M) by ag{v)w = gx{v,w) for v,w E T^M. This defines the gradient of a function / 
as grad/ = a~^{df). Define the divergence of a vector field X as the negative formal 
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adjoint of grad with respect to the scalar product {X,Y)x^2(^m) = Jm 9^^-^^'^) dvol(a;): 

(X,grad/)L2(jv/) = -(divX, /)i2(A^). 
The Laplace operator is then defined as (cf. [551156] ) 

A/ = divgrad/. (10) 

By definition this operator is self-adjoint and thus has only non-positive real eigenvalues 
by (A/,/)i2 = -(grad/, grad < 0. 

Based on the above approach, a corresponding energy form ("Dirichlet form", 
cf. [57]) can be defined 

£{u,v) = / (7a; (grad M, grad w) dvol(a;) = / gx{du,dv) dvo\{x), 
J M Jm 

which lends itself to a further way of defining a Laplace operator via the relation 

£{u,v) = -{Au,v)l2(^m)- (11) 

Geometrically, the Laplace operator measures the deviation of the function / from 
the mean value. More precisely, let S{x, r) = {y E M \ d{x, y) = r} denote the ball of 
radius r in M (in the Riemannian metric). Then 

where n denotes the dimension of the manifold M and a is the surface measure on 
S{x,r). This also motivates the definition of the Laplace operator as limit of finite 
difference operators 



A/(a;) = liml[ ^ w,f{p)-f{x)\ 



where Nr{x) is a finite set of points at distance r from x and Wp are suitably chosen 
weights. Such approximations to the Laplace operator are the basis of the method of 
finite differences in numerical mathematics. 

The Laplace operator can then be used to define a diffusion on M via the heat 
equation Au = dtu. The solution u{t,x) of the initial value problem ^(0,^) = f{x) 
defines a semi-group of operators Pj by 

u{t,x) = PJ{x). 

The semi-group property Ps+t = PsPt comes from the uniqueness of the solution u and 
translation invariance with respect to t of the heat equation. From the heat semi-group 
Pt the Laplace operator can be recovered as the infinitesimal generator 

A / = lim (13) 

which exists on a dense subspace of L^{M) under suitable continuity assumptions on 
the semi-group {Pt)t>o (cf. [58]). 
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In the fractal situation none of the above approaches can be used directly to define 
a Laplace operator. The main reason for this is that there is no natural definition of 
derivative on a fractal. But the above approaches to the Laplace operator on a manifold 
can be used in the opposite direction: 

• starting from a diffusion process that can be defined on fractals by approximating 
random walks. Then the Laplace operator can be defined as the infinitesimal 
generator. This is described in Sections 13.21 and 13.31 

• taking the limit of finite difference operators on graphs approximating the fractal 
gives a second possible approach to the Laplace operator, which is explained in 
Section [331 

• starting with a Dirichlet form S gives a third possible approach, which is presented 
in Section 1331 

3.2. Random walks on graphs and diffusion on fractals 

The first idea to define a diffusion on a fractal was to define a sequence of random 
walks on approximating graphs and to synchronise time so that the limiting process is 
non-constant and continuous. This was the first approach to the diffusion process on 
the Sierpihski gasket given in [HI [71 [8] and later generahsed to other "nested fractals" 
in [in]. Because of its importance for our exposition, we will explain it in some detail in 
this section. We will follow the lines of definition of self-similar graphs given in [591160] 
and adapt it for our purposes. 

We consider a graph G = (y{G), E{G)) with vertices V{G) and undirected edges 
E{G) denoted by {x, y}. We assume throughout that G does not contain multiple edges 
nor loops. For G C V{G) we call dG the vertex boundary, which is given by the set of 
vertices in V{G) \ G, which are adjacent to a vertex in G. For F C V{G) we define the 
reduced graph Gp by V{Gf) = F and {x, y} G E{Gf), if x and y are in the boundary 
of the same component of V{G) \ F. 

Definition 3 A connected infinite graph G is called self-similar with respect to F C 
V{G) and ^ : V{G) V{Gf), %f 

(i) no vertices in F are adjacent in G 

(a) the intersection of the boundaries of two different components ofV{G) \ F does not 
contain more than one point 

(Hi) (f is an isomorphism of G and Gp. 

A random walk on G is given by transition probabilities p{x, y), which are positive, if and 
only if {x,y} G E{G). For a trajectory (F„)„gNo of this random walk with Yq = x & F 
we define stopping times recursively by 

r„+i = min {k>Tm\Yk(^F\ {YtJ} , Tq = 0. 
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Then (lT™)mGNo is a random walk on Gp. Since the underlying graphs G and Gp are 
isomorphic, it is natural to require that {'^^^(yTm))mmo is the same stochastic process 
as (y„),igNo- This requires the validity of equations for the basic transition probabilities 

P {Yt„^, = ifiy) I Yt^ = ifix)) = P = y\Y„ = x)= pix, y). (14) 

These are usually non-linear rational equations for the transition probabilities p{x,y). 
The existence of solutions of these equations has been the subject of several 
investigations, and we refer to [6 n[62l[63[ [64] . 

The process (l^)neN() on G and its "shadow" (^„)neN() on Gp are equal, but they 
are on a different time scale. Every transition Yp^ — )■ on Gf comes from a path 

Yt„ — )■ Yt„+i ■■■—)■ yT„+i-i — ^ ^T„+i in a component of V{G) \ F. The time scaling 
factor between these processes is given by 

A = E(T„+i-T„)=E(ri). 

This factor is > 2 by assumption on F. More precisely, the relation between the 
transition time on Gp and the transition time on G is given by a super-critical (A > 1) 
branching process, which replaces an edge {(p{x), (p{y)} G Gp by a path in G connecting 
the points x and y without visiting a point in V{G) \ F (except for x and for y in the 
last step). 

In order to obtain a process on a fractal in R'', we assume further that G is 
embedded in M°' (i.e. V{G) C M°'). The self-similarity of the graph is carried over 
to the embedding by assuming that there exists a /3 > 1 (the space scaling factor) such 
that F = V{Gp) = f3V{G). The fractal limiting structure is then given by 

oo 

YG=[jp-^V{G). 

n=0 

Iterating this graph decimation we obtain a sequence of (isomorphic) graphs 
Gn = {f3^"'y{G),E{G)) on different scales. The random walks (Y'^^"^)^^^^ on G„ are 
connected by time scales with the scaling factor A. From the theory of branching 
processes (cf. p]) it follows that the time on level n scaled by A~" tends to a random 
variable. From this it follows that /^""^[a"] weakly tends to a (continuous time) 
stochastic process {Xt)t>o on the fractal Yq. 

Under the assumption that the graph G has an automorphism group acting doubly 
transitively on the points of the boundary of every component of G \ F (cf. [59|[60]). 
we consider the generating function of the transition probabilities Pn{x, y) = P(Fn = y \ 
Yq = x), the so called Green function of the random walk 

oo 

G{z \x,y) = ^pnix,y)z'^, 

n=0 

which can also be seen as the resolvent (/ — zP)^^ of the transition operator P. The 
replacement rules connecting the random walks on the graphs G and Gp result in a 
functional equation for the Green function 

G{z I ^{x), ^{y)) = f{z)Gi^{z) \x,y), (15) 
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where the rational function ipi^) is the probabihty generating function of the paths 
connecting two points vi,V2 G ip(y{G)) without reaching any point in ip{V{G)) \ {vi}. 
Then the time scahng factor is the expected number of steps needed for the paths 
counted hj tp, so we have 

A = V''(1). 

The rational function f{z) is the probability generating function of the paths starting 
and ending in vi G ip{V{G)) without reaching any other point in ip{V{G)). Equation 
(IT3]) becomes especially simple for a fixed point x of the map (p 

G{z \x,x) = f{z)G{jp{z) I x,x). 

It was proved in |60j that under our conditions on the set F the map can have at 
most one fixed point. 

The Koenigs function ^ oi ip around the fixed point z = 1 is given by 

^{\z) = ^($(z)), $(0) = 1, $'(0) = 1. 

This can be used to linearise this functional equation and to obtain precise analytic 
information about G{z \ x,x) (cf. [60l[65]). From this the asymptotic behaviour of 
the transition probabilities Pn{x,x) can be derived. In many examples these transition 
probabilities exhibit periodic fluctuations 

Pn{x,x) ~ n~— (cr(log;,n) + 0{n' )) , 

where a is a continuous, periodic, non-constant function of period 1 (cf. [60]). 

A first example of such graphs is the Sierpihski graph studied as an approximation 
to the fractal Sierpihski gasket. In this case we have for the probability generating 
function = ^ind A = 5. The random walk on this graph was studied in [8] in 
order to define a diffusion on a fractal set. Self similarity of the graph and the fractal 
have been exploited further, to give a more precise description of the random walk [66] 
and the diffusion [67]. In [60] a precise description of the class of graphs in terms of 
their symmetries is given, which allow a similar construction. In [68] this analysis was 
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carried further to obtain results for the transition probabihties under less symmetry 
assumptions using multivariate generating functions. 

In some examples (for instance the Sierpihski graph) it occurred that the function 
ip was conjugate to a polynomial p, i.e. 

" W 

which allowed a study of the properties of the random walk by referring to the classical 
Poincare equation. Properties of the Poincare equation have been used in [69] to study 
the analytic properties of the zeta function of the Laplace operator given by the diffusion 
on certain self-similar fractals. For more details we refer to Section 14.41 

Remark 3.1 The Sierpinski carpet is a typical example of an infinitely ramified fractal. 
In f7U^ approximations by graphs are used to define a diffusion on this fractal. By the 
infinite ramification, this approach is is more intricate than the procedure described here. 
By the very recent result on the uniqueness of Brownian motion on the Sierpinski carpet 
( cf. 138] ) this yields the same process as constructed by rescaling the classical Brownian 
motion restricted to finite approximations of the Sierpinski carpet in fW^lSTV 




Figure 4. The Sierpinski carpet 



3. 3. The Laplace operator as the infinitesimal generator of a diffusion 

Given a diffusion process {Xt)t>o on a fractal K we can now define a corresponding 
Laplace operator. At first we define a semi-group of operators At by 

AJ{x) = EJiXt) 
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for functions / G L'^{K). The semi-group property 

of the operators comes from the Markov property of the underlying stochastic process 
By [58| Chapter IX] this semi-group has an infinitesimal generator given as 

t^0+ t 

This limit exists on a dense subspace J-' of LP'{K) and is called Laplace operator on K. 
This name comes from the fact that for the usual Brownian motion on a manifold this 
procedure yields the classical Laplace-Beltrami operator. The function u{x,t) = Atf{x) 
satisfies the heat equation 

Au = dtU, u{x,0) = f{x). 

It was observed in the early beginnings of the development of the theory of diffusion 
of fractals that the domain of A does not contain the restriction of any non-constant 
differentiable function (cf. \7T\). 



3.4- The Laplace operator as limit of difference operators on graphs 

A totally different and more direct approach to the Laplace operator on self-similar 
fractals has been given by Kigami in [13]. The operator A is approximated by difference 
operators on the approximating graphs G*„. The graph Laplacians are given by 

^nf{x)= ^ p{x,y)f{y) - f{x) 

y~x,y£G„ 

as a weighted sum over the neighbours of x in Gn (y ~ x describes the neighbourhood 
relation in the graph In order to make this construction compliant with the 

approach via stochastic processes, these operators have to be rescaled appropriately. 
The correct rescaling is then given by the time scaling factor A introduced before, namely 

A/(x) = lim A"A„/(x). 

n— >oo 

3. 5. Laplace operators via Dirichlet forms 

Following the exposition in [I6l Chapters 2 and 3] we define a sequence of quadratic 
forms on the finite sets Vm given in Definition El 

Definition 4 Let V he a finite set. Then a bilinear form £ on liV), the real functions 
on V is called a Dirichlet form, if the following conditions hold: 

(i) Vu G liV) : £{u,u) > 

(a) £{u,u) = implies that u is constant on V 

(Hi) for u G (-iy) and u{x) = max(0, min(M(a;), 1)) the inequality £{u,u) > £(u,u) 
holds. 
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Definition 5 Let S be a Dirichlet form on the finite set V and let U he a proper subset 
ofV. Then the restriction of S to U is defined as 

'J^v,u{^){u, u) = mill {£{v, v) \ V E i(y),v\u = u} . (16) 

On the level of the coefficient matrices of the Dirichlet forms, the operation of restriction 
is given in terms of the Schur complement. 

Definition 6 Let {Vn, £n)n be a sequence of increasing finite sets K„ and Dirichlet forms 
En on Vn- The sequence is called compatible, if 

^y„+i,y„(^n+i) = En 
holds for all n. 

For a compatible sequence {yn,£n)n and a function / on K, the sequence 
(^n(/|y„, /|v„))n is increasing by definition and thus converges to a value in [0, oo]. 
This makes the following definition natural. 

Definition 7 Let {Vn,Sn)n be a compatible sequence of Dirichlet forms. Let 

V={f:K^R\\im£nif\v^J\vJ<oo]. 
Then for all f eV 

S{fJ)= Mm Sniflv.JlvJ 

defines a Dirichlet form on K, and T) is its domain. 

In order to make the sequence of Dirichlet forms coherent with the self-similar 
structure of K, we require the following self-similarity condition for £n 

m 

Sn+i{f, /) = A 5^ r-'Snif °F.-.fo F,), 

i=l 

where {i = 1, . . . , m) are positive weights and A is a proportionality factor. 
Furthermore, the sequence of forms has to be compatible, which amounts to the equation 

A7^y„yo o F,, ■ o F,))j =^o(-,-), (17) 

which comes as a solution of a non-linear eigenvalue equation. This equation plays the 
same role in the Dirichlet form approach as the equations ( I14p play for the approach 
via random walks. The Dirichlet form S on K is then defined as 

S{f, f) = lim A" Yl ^-'^o(/ oF^lv.Jo F^ 
where S" = {1, . . . , m}, = r^^ ■ ■ ■ for w = Wi . . . w„ and F^ = F^^ o ■ ■ ■ o F^^^ 



Remark 3.2 There are some additional technical problems concerning this construction 
of Dirichlet forms, which arise from the fact that in general the form is supported only 
on a proper subset of K . In [TSj Chapter 3] sufficient conditions for the weights Ti and 
the form Sq are given, which ensure that the form £ is supported on the whole set K . 
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Given a Dirichlet form £ together with its domain V and a measure /i on K, we 
can now define the associated Laplace operator on K by 

Wv eVnl^fi) : S{u,v) = -{A^u,v)mf,), (18) 

which defines A^, if V is dense in L'^{fJ.). Notice that this is the same equation as in the 
manifold case (fTTI) . 

In the case of a self-similar fractal K as described in Section[2]the "natural" measure 
on K is the according Hausdorff measure H^. In this case we omit the subscript /i. 

For more information on Dirichlet forms in general and their applications to the 
description of diffusion processes, we refer to the monograph [57]. For the specifics of 
Dirichlet forms on fractals we refer to [T61I72] . 



4. Spectral analysis on fractals 
4-1- Spectral analysis on manifolds 

Let us start with a short discussion of the manifold case, which is again somehow 
complementary to the fractal case. As described in Section [3] the Laplace operator on 
a manifold M is defined via a Riemannian metric g. If M is compact with smooth (or 
empty) boundary (for simplicity) the Laplace operator has pure point spectrum. The 
eigenvalues — are all real (by self-adjointness) and non-positive (by definition (fTOj) ). 
Denote the normalised eigenf unctions by ipk- Then the heat kernel can be written as 

oo 

pt{x,y) = Y,e-^'''M^)My)- (19) 

fc=0 

This expression yields the trace of the heat kernel 

» oo 

K{t)= / pt(x,x)dvol(x) = Ve-^^*. (20) 

On the other hand, if M is closed, the asymptotic behaviour of the heat kernel for 
t — )■ 0+ can be described very precisely 



oo 

p,{x,y) = (47rt)-te-^(-'^)'/(^*) ^t^a,(x,i/), (21) 

£=0 

where the functions ae can be computed iteratively by solving certain second order 
partial differential equations (cf. [73l[71]). In particular, ao{x,y) = 1. The values 
ae{x,x) can be expressed in terms of the functions Ay d{x,y)'^"^\^_^ (cf. [75]). In the 
case of manifolds M with smooth boundary, further terms involving half integer powers 
of t occur 

oo oo 

Pt{x, x) = t~^ ^ t^ai{x, x) + t~^ ^ t2bi{x, x); (22) 

£=0 i=l 

the terms encode curvature information of the boundary dM. For more details we 
refer to 
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Especially, this gives an asymptotic expansion of K{t) (in the case of a closed 
manifold) 

K{t) = {ATTty'^^t^ ae{x,x)dYo\{x). (23) 

From the first order asymptotic relation K{t) ~ vol(M)/(47rt)t the asymptotic 
behaviour of the counting function 

Nix) = 1 (24) 

\k<x 

can be obtained by a Tauberian argument giving the classical Weyl asymptotic relation 

, VOl(M) n . n 

N{x) = — ^-^x^ + o(x2 ), (25) 

where Bn denotes the volume of the n-dimensional unit ball. The eigenvalues of A 
constitute the frequencies of the oscillations in the solutions of the wave equation 
Au = Utt- This led to M. Kac's famous question "Can one hear the shape of a drum?" 

(cf. m)- 

n — l 

The precise order of magnitude of the error term has been determined as x^~ in 
the case of a smooth boundary of M in [78]. In the case of fractal boundary, this is 
subject of the Weyl-Berry conjecture stated in [79|l80] . The original conjecture estimated 
the error term by 0{x 2 ), where dnidM) denotes the Hausdorff dimension of the 
boundary of M. This was shown to be false in [81]. The Hausdorff dimension was then 
replaced by the Minkowski dimension dM{dM). It was shown in [82] that the error term 
is 0{x2) for any d > dM{dM). Furthermore, it was proved in [8^ that the modified 
conjecture is true for dimension n = 1; in [HI] counterexamples for dimension n > 1 are 
presented. In these counterexamples disconnected manifolds are considered, thus the 
conjecture could still hold for connected manifolds M with fractal boundary. 

More precise information on the eigenvalues is contained in the spectral zeta 
function 

C^(,)= ^(A,)-, (26) 

the Dirichlet generating function of {Xk)k- The zeta function is connected to the trace 
of the heat-kernel by a Mellin transform 

CA(s)r(s) = / K{t)t'-^dt, for3?s>-. (27) 
Jo 2 

The right hand side has an analytic continuation to the whole complex plane, which can 

be found by using the asymptotic expansion ( 123|) . Furthermore, (in the case of a closed 

manifold) the right hand side has at worst simple poles at the points s = — i (with 

£ G {0, 1, . . . , [^^J for even n and £ G Nq for odd n) with residues 



Res Uis) = ae{x,x)dvol{x) 
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and special values (for £ G N) 

Ca(-^) = 



for n odd 

{—lyU Jj^ ae{x, x) dvo\{x) for n even 



and 



Ca(0) 



— dim ker A for n odd 

/jv/ {x, x) d vol(x) — dim ker A for n even 

(cf. L73J). The special value Ca(0) '^^^ be interpreted as the negative logarithm of 
the determinant of A. In the case of a manifold M with smooth boundary a similar 
reasoning gives simple poles of the zeta function at the negative half-integers — i, — |, . . . 
and at the points | — — |,..., '^^^^^^ with residues depending on the functions bi. 

In the case of a fractal, the above approach to the spectral zeta function, its analytic 
continuation, and its finer properties can not be used. The reason for this is that no 
asymptotic expansion of the heat kernel is known in the fractal case. In Section we 
comment on the known upper and lower estimates for the heat kernel. 

On the other hand, for fractals having spectral decimation, the eigenvalues can be 
described very precisely. In this case they turn out to constitute a finite union of level 
sets of solutions of the classical Poincare functional equation. This allows then to obtain 
the analytic continuation of the spectral zeta function from the precise knowledge of the 
asymptotic behaviour of the Poincare function; the asymptotic expansion of the trace of 
the heat kernel can then be obtained from the poles of the zeta function, reversing the 
argument in (!27|) . The eigenvalue counting function can then be related to the harmonic 
measure on the Julia set of a the polynomial governing the spectral decimation. 



4-2. Spectral decimation 

It has been first observed by Fukushima and Shima [T71|20l[2T] that the eigenvalues 
of the Laplacian on the Sierpihski gasket and its higher dimensional analogues exhibit 
the phenomenon of spectral decimation. Lateron, spectral decimation for more general 
fractals has been studied by Malozemov, Strichartz, and Teplyaev [2211231185] . 

Definition 8 (Spectral decimation) The Laplace operator on a p. c. f. self-similar 
fractal G admits spectral decimation, if there exists a rational function R, a finite set 
A and a constant A > 1 such that all eigenvalues of A can be written in the form 

lim A"i?(-")({w}), weA.men (28) 

n— >oo 

where the preimages of w under n-fold iteration of R have to be chosen such that the 
limit exists. Furthermore, the multiplicities Pm{w) of the eigenvalues depend only on w 
and m, and the generating functions of the multiplicities are rational. 

The fact that all eigenvalues of A are negative real implies that the Julia set of R has 
to be contained in the negative real axis. We will exploit this fact later. 
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In many cases such as the higher dimensional Sierpihski gaskets, the rational 
function R is conjugate to a polynomial. The method for meromorphic continuation 
of (a given in Section 14.41 makes use of this assumption. Recentlty, Teplyaev [85] 
showed under the same assumption that the zeta function of the Laplacian admits 
a meromorphic continuation to 3fts > —e for some e > depending on properties of the 
Julia set of the polynomial given by spectral decimation. His method uses ideas similar 
to those used in |67] for the meromorphic continuation of a Dirichlet series attached to 
a polynomial. Complementary to the ideas used here, Teplyaev's method carries over 
to rational functions R. 

4-3. Eigenvalue counting 

As in the euclidean case the eigenvalue counting function 
Nix) = J2 1 

Xk<x 

measures the number of eigenvalues less than x. In a physical context this quantity 
is referred to as the "integrated density of states". It turns out that N{x) in many 
cases does not exhibit a pure power law as in the euclidean case, but shows periodic 
fluctuations. One source of this periodicity phenomenon is actually spectral decimation, 
especially the high multiplicities of eigenvalues, as will become clear in Section 14.41 

Recall the definition of harmonic measure on the Julia set of a polynomial p of 
degree d (cf. [86]): the sequence of measures 

^" = ^ Yl 

p(")(x)=5 

converges weakly to a limiting measure /i, the harmonic measure on the Julia set of p. 
The point ^ can be chosen arbitrarily. 

Assume now that the Laplace operator on the fractal K admits spectral decimation 
with a polynomial p. Then the relation lim„^oo X'^p^~'^\{w}) G -B(0, x) can be translated 
into 

lim p^"'\X~"z) = w and \z\ < x. 

n— >oo 

Here B{0,x) denotes the ball of radius x around 0. By general facts about polynomial 
iteration (cf. [26]), the limit exists and defines an entire function of z, the Poincare 
function ^{z). The number of eigenvalues with m = in B{0,x) is the equal to 

n{x) = j2nm 

with 

N^{x) = 4^{zeB{0,x) I <!>{z)=w}. 
The following relation can be obtained from the definition of harmonic measure 
lim (i-"Ar^(A"x) = /i ($-^(5(0, x))) ; 
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this holds for all x small enough to ensure the existence of the inverse function $ ^ on 
5(0, x). 

In [871 Theorem 5.2] we could prove a relation between the asymptotic behaviour 
of the partial counting functions N^{x) and the harmonic measure of balls fi{B{0,x)). 
The existence of the two limits (p = log^ d) 

lim x~^Nw{x) 

x—>oo 

^mt-ffi{B{0,t)) 

is equivalent. We conjectured there, that these limits can only exist, if p is either a 
Chebyshev polynomial or a monomial. These are the only cases of polynomials with 
smooth Julia sets (cf. [88]). 

Summing up the above discussion, the eigenvalue counting function can be written 

as 

oo 

^(^) = Y.Y1 PmHNU^'"'^)- (29) 

Notice, that for fixed x, these sums are actually finite. In the known cases, such as 
the Sierpihski gasket, the growth of (3m{w) is stronger than d"^, which implies that the 
terms for large m (with Nuj{X~"^x) still positive) become dominant in this sum. This 
shows that multiplicity of the eigenvalues has the main influence on the asymptotic 
behaviour of N{x). Furthermore, this explains the presence of an oscillating factor in 
the asymptotic main term of N{x). We will discuss that in more detail in Section 



4-4- Spectral zeta functions 

As in the euclidean case, the eigenvalues of the Laplace operator A can be put into 
a Dirichlet generating function. This will later allow to use methods and ideas from 
analytic number theory to obtain more precise asymptotic information on N{x). The 
zeta function is again given by 

Ca(^) = Yl 

where all eigenvalues are counted with their multiplicity. The zeta function is related 
to the eigenvalue counting function by 

/•oo /"OO 

Ca{s) = / x-'dN{x) = s N{x)x-'-^dx. 
Jo Jo 

The second relation identifies (a as the Mellin transform of the counting function N{x). 

In the sequel we will exploit the consequences of spectral decimation. Not too 
surprisingly after Definition [8], iteration of polynomials will play an important role in 
this discussion. Furthermore, since the relation ( 128|) can be expressed in terms of the 
Poincare function /, properties of this function will be used to derive the meromorphic 
continuation of Ca to the whole complex plane. 
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Under the assumptions of spectral decimation, the Juha set of the polynomial p is 
a subset of the non-positive reals, which contains 0. By [87, Theorem 4.1] this implies 
that A = p'{0) < <f. By [871 Theorem 4.1] equality can only occur, if p is a Chebyshev 
polynomial, which would correspond to spectral decimation on the unit interval (viewed 
as a self-similar fractal). Thus in the cases of interest, we have that p = logx d < ^- The 
Poincare function $ is then an entire function of order p. 

In order to find the analytic continuation of Ca(s) to the whole complex plane, we 
analyse the partial zeta functions 

a,«.(5)= J2 (30) 

Since = 1 — ^"^^ is a function of order p = log^^ci < it can be expressed as a 
Hadamard product 



for w; = we have the slightly modified expression 



${-m)=0,/.7^0 

Taking the Mellin-transform of log^w, which exists for — 1 < 3fJs < —p, we get 

poo ^ 

M^(s)= / (log$^(x))x^-Mx = — (<s>A-s). (31) 

Thus for finding the analytic continuation of C<!>^w{s) to the left of its abscissa of 
convergence, it suffices to find the analytic continuation of My^^s) for > —p. 
Following slightly different lines as in [69], we consider the function 

pmz))-w ^ ^^jxz) 

Then we have 

log^t„(z) = \og^^{Xz) -dlog$^(z) - log ad - {d- 1) log(-w) 

and this function tends to exponentially for z +oo. Taking the Mellin transform 
we obtain 

{X'-d)MUs) 

POO 

= / (log<l>^(Ax) - rflog$u,(x)) x^'^Mx 
Jo 

for - 1 < 3fJs < -p 

= / {log ^ u]{Xx) — d log ^y^{x) —log ad — {d — l)\og{—w)) x^~^ dx 
Jo 

-I- (logad+ (rf- l)log(-u;))- + / (log<l>^(Ax) - rflog$^(a;))a;"~Mx 

s Ji 

for 3fJs > -1 
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(log ^^,{x)) x^''^ dx 

/oo 
(log$^(Ax) -c^log$^„(x) -logorf- {d- l)\og{-w))x'~^ dx 



for 3fJs > 

(log\I>ti,(x)) x"""-^ dx. 



Reading the fourth hne of this computation shows that M^(s) has a simple pole at s = 
with residue 

ResM^(s) = i^ + log(-w). 

Furthermore, this computation shows that M^{s) is holomorphic in the half-plane 
3fJs > 0. Using fl3T]) gives the analytic continuation of Cq,^w{s) for JJs < 

C^As) = r {log^M))x~'-'dx. (32) 



7r(l - dX\ _ 

This shows that C'S>,w{—^) = for m G Nq (for s = the double zero of ssinvrs cancels 
the simple pole of M^(— s)). These could be called the "trivial zeros" as in the case of 
the Riemann zeta function. Furthermore, we obtain 

C$,«,(o) = ^— Y + iog(-w). 

The equation (152]) even lends itself to the numerical computation of values of C$,«)(s) 
for 3fJs < 0, as we will se in Section I4.5.2I 

By our assumption on spectral decimation, the generating functions of the 
multiplicities of the eigenvalues are rational 



m=0 

Thus we can write the spectral zeta function of A as 



oo 



Ca(s) = EE/^-H E (A»'^ = $^i?.(A-^)a,.(.). (33) 

Since all functions involved in the last (finite) sum are meromorphic in the whole complex 
plane, we have found the meromorphic continuation of Ca to the whole complex plane. 
The functions C'S>,wis) have only simple poles in the points s = log^^ + |^ (k G Z). 
Furthermore, the residues of these poles do not depend on w. All other poles of Ca come 
from the poles of the functions -Ru;(A~*). Since these are rational functions of A"**, their 
poles are equally spaced on vertical lines. 

Summing up, Ca('S) has a meromorphic continuation to the whole complex plane 
with poles in the points 

2kTTi 

s = - log;^ Pujj + T With A; G Z, 

log A 
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where the are the poles of the rational functions Rw{z) and (at most) simple poles 
in the points 

s = log^ d + with k eIj. 

log A 

Furthermore, since the functions C'^^w{s) are bounded for 3fts > log;^ c? + e and the 
functions Ryj{\~'^) are bounded along every vertical line, which contains no poles, the 
function Ca('S) is bounded along every vertical line c + it for c > log;^ d, which does not 
contain a pole of any of the functions Ryj{\~'^). 

In the case of the Sierpihski gasket and its higher dimensional analogues, the 
rightmost poles of Ca come from the rational functions i?^(A~'^). Furthermore, the 
relation 

^i?^(l/rf) = 

holds, which amounts in a (still mysterious) cancellation of the poles of the functions 
C*,™ in fl33l) . Thus, the analytic behaviour of the function is mainly governed by the 
functions Ryj. This means that in this respect the effects of multiplicity prevail over the 
individual eigenvalues. 

We now investigate the asymptotic behaviour of N{x) under the assumption that 
Pm{w) grows exponentially faster than d"^ for some w & A. This implies that the 
corresponding term i?s(A~*) has poles to the right of 3fJs = log;^ d. We use the classical 
Mellin- Perron formula (cf. [89j) to express N{x) in terms of Ca{s) 

1 r-^'^ ds 

for any c such that -R^(A~'*) has no poles in the half-plane 3fJ(s) > c. Now the line of 
integration can be shifted to the left to 3fJ(s) = c' for A'^' > d but such that at least one 
of the functions -R^(A~'^) has poles to the right of c'. This is justified, because Ca(<7 + ^^) 
remains bounded for |t| — )■ oo and a > c' . Then we have 

(c'+iT 
— I Ca(s)x^-- y Res 
/_iT 27r|/c|<TlogA 

where we denote the real part of the poles of i?^(A^*) by (is/2 (the "spectral dimension"). 
Now the limit of the integral J^,^-)'" can be shown to exist for T — > oo, which shows that 
the limit 

hm y Res 

T^oo ^-^ s=ds/2+2kni/logX S 
27r|fc|<riogA 

exists. This can be rewritten '^s/2i7(log^ x) for a periodic continuous function H 
given by its Fourier expansion. The limit of the integral can be shown to be 0{x'^'). 
Thus we have shown 

N{x) = x'^s/2^(iog^ x) + Oix"'). (34) 
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Especially, the limit \imx^^x~'^^^'^N{x) does not exist. 

We remark here that there exists totally different approach to zeta functions of 
fractals due to M. Lapidus and his collaborators [2HI EHl EDI EI] . In this geometric 
approach the volume of tubular neighbourhoods of the fractal G 



is analysed. The asymptotic behaviour of Vg{b) for e — )■ gives rise to the definition of 
a zeta function. In this geometric context the complex solutions of (jl]) occur as poles 
of the zeta function; they are called the "complex dimensions" of the fractal G in this 
context. This approach is motivated by the definition of Minkowski content, which itself 
turned out to be too restrictive to measure (most of the) self-similar fractals. 

4-5. Trace of the heat kernel 

The diffusion semi group At introduced in Section 13.31 can be given in terms of the heat 
kernel pt{x,y). As opposed to the situation in the euclidean case, the knowledge on 
behaviour of the heat kernel for t — )■ is by far not as precise. Kumagai ^92j proved the 
following lower and upper estimates of the form 



where ds and c?^ are the spectral and the walk dimension of the fractal, respectively. 
These dimensions are related to the Hausdorff dimension df of the fractal via the so 
called Einstein relation dgd^ = 2df. In the fractal case usually d^ > 2, as opposed to 
the euclidean case, where d^ = 2, which implies ds = df. 

Only recently, a conjecture by Barlow and Perkins [8] could be proved by Kajino 
[23] • Namely that for any x E K the limit 



does not exist for a large class of self-similar fractals. This gives an indication, why 
obtaining more precise information on the heat kernel than the estimates (1351) and fl36l) 
would be very difficult. 

Even if the precise behaviour of the heat kernel seems to be far out of reach, the 
trace of the heat kernel 



J K 

can still be analysed in some detail under the assumption of spectral decimation. The 
reason for this are the two relations between K{t), N{x), and Ca 



Vcie) = vol {{x e R" I d{x, G) < e}) 




(35) 



(36) 



limt 2 pt(x, x) 





(37) 



(38) 



Laplace Operators on Fractals 



26 



The first expresses K{t) as the Laplace transform of N{x), the second gives Ca as Mellin 
transform of K{t). 

Given the precise knowledge on the zeta function obtained in Section 14.41 we can 
use the Mellin inversion formula to derive asymptotic information on K{t) for t — )■ 0. 
We have for c > ^ 



where integration is along the vertical line = c. We notice that the Gamma 
function decays exponentially along vertical lines, whereas Dirichlet series grow at most 
polynomially by the general theory of Dirichlet series (cf. [Hj)- Thus convergence of the 
integral is guaranteed. 

Shifting the line of integration in f l39p to 3?s = — M — | (Af 6 N) and taking the 
poles of the integrand into account, we obtain for t — t- 



Here, H and Hj are periodic continuous functions of period 1, whose Fourier coefficients 
are given as the residues of (/:^{s)T{s) in the poles on the lines 3fJs = ^ = % 

respectively. By the strong decay of the Gamma function, these functions are even real 
analytic. The values aj + 2km/ log A come as the poles of the functions Ryj{\~'^). This 
shows that there are only finitely many aj. Since M can be made arbitrarily large, the 
error term decays faster than any positive power of t for t — )■ 0. 

The existence of complex poles of the zeta function thus implies the presence of 
periodically oscillating terms in the asymptotic behaviour of the trace of the heat kernel 
for t — )■ 0. This implies that the limit Yim.t^Qt'^^/'-^pt^x^x) does not exist on a set of 
positive measure for x in accordance with the above mentioned result by N. Kajino [93] . 
The consequences of the presence of complex poles of the zeta function to properties 
of the heat kernel, the density of states, and the partition function as well as the 
physical implications of the resulting fluctuating behaviour in these quantities have 
been discussed in [30|[3T]. 

4-5.1. Casimir energy on fractals As an application of the spectral zeta function and 
its properties we compute the Casimir energy of a fractal. We follow the lines of the 
exposition in [76] and refer to this book for further details. 

Consider the differential operator P = — ^ + A on (K/^Z) x G. As usual the 
parameter /3 = 1/kT. The zeta function of the operator P is then given by 




(39) 



C— lOO 




(40) 




Using the theta function relation 
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we obtain 

The derivative Cp(0) can be interpreted as — logdet P, the logarithm of the (regularised) 
determinant of P. 

For the rescaled operator P/ /i^ we have 

Cp/,2(0)=CM0) + Cp(0)ln/x2. 

From fHTl) we compute 

Cp(0) = -/3 Res Ca{s), 

which vanishes in the case of the Sierpihski gasket, because Ca does not have a pole at 
— ^. Furthermore, we get 

The integral and the summation over n can be evaluated explicitly, which finally gives 

00 



(1 \ 
^ .7=1 



e 

The case of the Sierpihski gasket and other self-similar fractal shows an important 
contrast to the manifold case, where the zeta function (generically) has a pole at — | 
(see Section l^?T|) ■ In the manifold case the value Ca(— |) has to be replaced by the finite 
part of (/\{s) at s = the function minus the principal part at the polar singularity 
(cf. [76]). 

Now the energy of the system is given by 



-[ eV^. - 1 

Letting /3 — )■ 00, which is equivalent to letting temperature tend to 0, gives the Casimir 
energy 

Ecas = -^Ca (^~2 
for fractals, whose zeta function has no pole at — ^. 



4-5.2. Numerical computations We will now describe the numerical computation of 
the values C$,«;(— 1/2), which are needed for the computation of Ca(— 1/2) in the case 
of the Sierpihski gasket. This fractal admits spectral decimation with the polynomial 
p{x) = x{x + 5) (cf. [SlITT]). The corresponding Poincare function is then given by the 
unique holomorphic solution of the equation 



$(5^) = $(z)($(z) +5), <l'(0) = 0, $'(0) = 1. 



(42) 
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We use the expression for Ca(s) for Dirichlet boundary conditions derived in \W\ 
Section 7] 

C2is) = + 5.(5._')(5._3) C^.-3(^) + (5^ -1x5" -3) ^^-"^^^^- 

Combining this with ( 132|) and inserting s = — ^ gives 

155 + 135^5 p, ^ 245 + 47^5 ^ / x -i ^ 1 

H J logW„3(x)x 2 dx — J logW_5(x)x 2 dx . 

From the general information on the asymptotic behaviour of Poincare functions 
derived in |87l|95] we obtain estimates of the form 



exp(Cix'') < < exp(C2a;0 

for positive constants Ci and C2, p = logs 2 valid for x > Xq. Such estimates can 
be proved easily by first showing the estimate for an interval of the form [xq, 5xo] and 
then extending it by using the functional equation fH2l) . For instance, we used Ci = 1, 
C2 = 1.08, and Xq = 10. Then the functions log\l/^(x) tend to like exp(— C2X^) for 
X -7- 00. Given a precision goal e > 0, we choose T so large that 



/oo ^ 
exp(— C2X^)x^2 dx < e. 



Then the improper integrals in fH3|) can be replaced by . The functions \l/«,(x) can be 
computed to high precision by using the power series representation of $(x) for |x| < 1 
and the functional equation fj42|) to obtain 

=p(^'+i)($(x/5'=+^)) 

for 5*^ < |x| < 5^"^^. This allows the computation of the remaining finite integrals up to 
precision e. We obtained 

E^^ = 0.5474693544... 

for the Casimir energy of the two-dimensional Sierpihski gasket with Dirichlet boundary 
conditions. 

Similarly, we have for Neumann boundary conditions (cf. [69]) 

^a(^) = _ ^1(3. _ 3) ((2 ■ 5^ - 5)a,-3(s) + 5^a,-5(s)) . 

This gives by the same numerical estimates as before 

E^a, = 2.134394089264... 

for the Casimir energy two-dimensional Sierpihski gasket with Neumann boundary 
conditions. 
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4-6. Self- similarity and the renewal equation 

Let K he a post-critically finite self similar fractal with self-similar structure 
{K, S, (-Fj)™ ^) with contraction ratios and Hausdorff-dimension p. Let K be equipped 
with a Dirichlet form with parameters rj and A as described in Section 13.51 Denote by 
N-d{x) and Ntss{x) the eigenvalue counting functions for the Laplace operator under 
Dirichlet and Neumann boundary conditions, respectively. Then the following crucial 
fact was observed in [96] 

Y^N^[p^x\ <N^{x)<N^{x)<Y,N^\p^xY (44) 

where ai are the contraction ratios introduced in Section [2], p is the Hausdorff dimension 
of and the are the weights used for the construction of the Dirichlet form E in 
Section [3751 Furthermore, the inequalities 

iVD(2;)<iVN(x)<iVD(x) + #(K,) 

hold, where Vq is defined in Definition |2l Setting 



p\ 1/2 



7. . ^ 
we end up with the equation 

iVD(x) = ^iVD(7^)+^?(x), (45) 

i=l 

where g[x) is defined as the difference between the left hand side and the sum on the 
right hand side, and remains bounded by ( H^ . A similar equation holds for A^N(a;). 

Equation (H5i) can be transformed into the classical renewal equation occurring in 
probability theory (cf. [97]). The asymptotic behaviour of the solutions of this equation 
is described by the following theorem (stated as in [T6]). 

Theorem 4.1 ("Renewal Theorem") Lett* > and f be a measurable function on 
M such that f{t) = for t < t* . If f satisfies the renewal equation 

N 

where ai, . . . , Oat are positive real numbers, pj > for j = I, . . . , N and Yl'j^=iPj — ^■ 
Assume that u is non-negative and directly Riemann integrable on M with u{t) = for 
t < t* . Then the following conclusions hold: 

(i) arithmetic (or lattice) case: if the group generated by the aj is discrete (i.e. there 
exist aT > and integers rrij with greatest common divisor 1 such that aj = Trnj; 
all the ratios ai/aj are then rational), then limj^oo |/(^) ~ G{t)\ = 0, where the 
T -periodic function G is given by 

/ N \ ^1 oo 
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(a) non- arithmetic (non-lattice) case: if the group generated by the aj is dense in M 
(i.e. at least one of the ratios Oi/aj is irrational) , then 



N ^ -1 



}^f{t)= [f^P.aA j u{t)dt 



Corollary 4.1 Let f he a solution of the equation 

m 

1=1 

where < 7^ < 1 and g is a bounded function. Let ds be the unique positive solution of 
the equation 

m 

i=l 

then the following assertions hold: 

(i) arithmetic (or lattice) case: if the group generated by the values log7j is discrete, 
generated by T > 0, then 

f{x)=x''^/^ (G((logx)/2) + o(l)) 

for a periodic function G of period T . 

(a) non- arithmetic (or non-lattice) case: if the group generated by the values log'jj is 
dense, then 

lim /(x)x-'^s/2 

exists. 

The corollary is an immediate consequence of the Theorem by setting /(e*) = e'^s*/^-F(t) 
and applying the Theorem to F and g{e*)e~^^*^'^. 
Summing up, we have the following theorem. 

Theorem 4.2 ( |96L Theorem 2.4]) Let (i^', S, (Fj)™ be a self similar structure 
with contraction ratios Ui and Haus dor ff- dimension p. Let K be equipped with a Dirichlet 
form with parameters ri and X as described in Section lSTSi Let ds be the unique positive 
solution of the equation 

i=l ^ ' 

the "spectral dimension" of the harmonic structure of K. Then the following assertions 
hold: 

(i) lattice case: if the group generated by the values log(rja;f/A) is dicrete, then 
Nuix) = x"^'^ (G((log(x))/2) + 0(1)) 
for a periodic function G, which is non-constant in general. 
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(a) non-lattice case: if the group generated by the values \og{ria^ / X) is dense, then 
lim A/'D(x)a;-'^s/2 

x—^oo 

exists. 

The behaviour of N^{x) for x oo is the same. 

Remark 4.1 Similar ideas are used in 19^ to study the asymptotic expansion of the 
eigenvalue counting function of the Laplace operator on an open set (3, which is formed 
from an open set Gq with smooth boundary as 

oo 

e = u ^« 

n=0 

with 

m 

Gn+1 = 

i=l 

with similitudes Fi and all the unions assumed to be disjoint. 

Remark 4.2 Very recently, N. Kajino f9^ could extend this approach to non-p. c. f. 
fractals such as the Sierpinski carpets. He observed that in order to have an inequality 
of the form ( f^^ it suffices to have the corresponding self- similarity property 

m 

S{f,f) = XY,rr'S{foF,,foF,) 

i=l 

of the underlying Dirichlet form. This self- similarity together with symmetries 
characterises the Dirichlet form on the Sierpinski carpet uniquely, as was shown in (3^ 



5. Exploiting self-similarity 

5.1. Decimation invariance and equations with rescaling 

As has been described in Section [3121 diffusion processes on decimation invariant finitely 
ramified fractals, belonging to a certain class, may be defined as limits of discrete simple 
symmetric nearest-neighbour random walks on the associated approximating graphs. 
The analysis is based on a renormalisation type argument, involving self-similarity and 
decimation invariance . 

Consider, for example, the Sierpinski gasket F. It can be approximated by 
"Sierpinski lattice" graphs Let X^^"^ be a simple random walk on G„. Then, 

according to Goldstein [6] and Kusuoka [7j, 

as n — )■ cxD, where Xt is a diffusion process on F. ( Here, =^ means convergence in 
distribution and [.] is the integer part function.) 

Xt is a Markov process with continuous sample paths (in fact, a Feller process), 
which is itself invariant under the rescaling x — )■ 2x, t — 5t. 
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As has been explained in Section I3.2[ equation f llSp in this instance has the form 
<l>(Az) = i?($(z)), (46) 

2 

where A = 5; R{z) = ipi^z) = -^-^ is a rational function of z. 
Denoting here ^{z) = we have also 

^(A;z) = P{^{z)) (47) 

where A = 5; P{z) = Az"^ — 3^; is a polynomial. 

Both fj^6|) and fj47|) are examples of the Poincare equation which will be discussed 
below. 



5.2. Functional equations in the theory of branching processes 

Iterative functional equations, and the Poincare equation, in particular, occur also in 
the context of branching processes (cf. [9]). 
Here a probability generating function 

oo 
n=0 

encodes the offspring distribution, where Pn > is the probability that an individual 
has n offsprings in the next generation (note that g(l) = 1). The growth rate A = q'{l) 
determines whether the population is increasing (A > 1) or dying out (A < 1). In 
the first case the branching process is called super- critical. The probability generating 
function q^"\z) (n-th iterate of q) encodes the distribution of the size X„ of the n-th 
generation under the offspring distribution q. In the case of a super-critical branching 
process it is known that the random variables A~"X„ tend to a limiting random variable 
Xoo- The moment generating function of this random variable 

f{z) = Ee~'^^ 

satisfies the functional equation (cf. [9]) 

f{Xz) = q{f{z)), (48) 

which is yet another example of the Poincare equation. 

If q{z) = p(l/^-^ , where P is polynomial, then f HSj) coincides with P7|) . 

In general, such a reduction to P7|) is possible, when q is conjugate to a polynomial 
by a Mobius transformation. 



5.3. Historical remarks on the Poincare equation 

In his seminal papers |100lll01j . H. Poincare has studied the equation 

f{Xz) = R{f{z)), zee, (49) 

where R{z) is a rational function and A G C. He proved that, if -R(O) = 0, -R'(O) = A, and 
|A| > 1, then there exists a meromorphic or entire solution of (l49l) . After Poincare, (l49l) 
is called the Poincare equation and solutions of (H9l) are called the Poincare functions 
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. The next important step was made by Valiron [102(1103] . who investigated the case, 
where R{z) = P{z) is a polynomial, i.e. 

fiXz) = P{fiz)), zee, (50) 

and obtained conditions for the existence of an entire solution f{z). Furthermore, he 
derived the following asymptotic formula for M(r) = ma:X\z\<r \f{z)\'- 

In M(r) -r'^g (-^) , r ^ oo. (51) 



InlAI 



Inm 



Here Q{z) is a 1-periodic function bounded between two positive constants, p — 
and m = deg P{z). 

Various aspects of the Poincare functions have been studied in the papers 
[871I951 ITMIT051IT061I1071I1081IT09] . 



5.4- Applications: diffusion on fractals 

In addition to ( ISTll . in applications (see Sections 14.41 and 14.51 above, in particular ) it is 
important to know asymptotics of entire solutions f{z) in certain angular regions and 
even on specific rays re^^ of the complex plane. 

In the following section we present some recent results in this direction ( [STflOS]). 

Note, that in the aforementioned applications ( diffusion on fractals ) scaling factor 
A is real and, also, polynomial P{z) has real coefficients. 

Nevertheless, we prefer to start our exposition in section 1531 below from the general 
case and turn to the real case afterwards. We think that, these general facts on the 
Poincare equation are interesting in their own, and hope that, maybe, they will also find 
applications in the future. 



5.5. Asymptotics along spirals and asymptotics along rays 

Here we intend to describe further results of Valiron's type . We start our presentation 
from some of our recent results [871195]. 

It turns out that asymptotic behaviour of the Poincare functions heavily depends 
on the arithmetic nature of argA, where A is the scaling factor in (!50l) . Denote 
arg A = 271 (3 . 

The following statement is true, for irrational /3 ( [95]) : 

Theorem 5.1 //argA = 27r/3 and (5 is irrational, then f{z) is unbounded along any ray 
d. Moreover, if we denote f{z) = ln\f{z)\ (where main branch of logarithm is taken) 
then there exists a sequence r„ — )■ oo, such that the limit 

hm = L (52) 

n—>-oo Tn 

exists and L > 0. 
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As far as we know this phenomenon has not been mentioned in the hterature before. 

On the other hand, if /3 is rational (and, in particular, if /3 = 0, i.e. A is real) f{z) 
may be bounded on some rays and even in whole sectors. Nevertheless, for rational 
(3, the limit ( 152|) still exists under some additional assumptions. Denote (3 = t/s and 
suppose that t, s are relatively prime. Put q = (note that 1 < g e R). 

The following result is true for rational /3 and, in particular, for real A ( |95j): 

Theorem 5.2 Suppose that either |A| > or s > 2p. Then f{z) is unbounded on any 
ray, and one can find a geometric progression Tn = Q^tq, [vq > 0), for which the limit 
[5^] exists and L > 0. 



The above results are based on the following theorem on asymptotics of f{z) along 
spirals (geometric progressions) of the form Zn = X'^zq, where zq = roe*^" G C will be 
specified below. 

Consider the Poincare equation 

f{Xz) = P{f{z)), (53) 



where 



P{Z) = Z^+ Pm-lZ""-' + ...+P1Z+P0. (54) 



Denote 



K = max{\po\,. . .\p^_i\}. (55) 
Suppose that zq = roc*^" G C is such a point, that 

\f{zo)\>max{e,2mK}. (56) 
Let us use the following notations 

z„ = A%; ^(;^) = log|/(z)|; P = (57) 

log |A| 

in part already introduced earlier. (It follows from fl56|) that (p{zq) > 1 ) 

Theorem 5.3 Suppose that / f55|) is satisfied. Then the limit along the spiral Zn = X'^zq 

hm ^ = L{zo) (58) 



exists and 



^{z,)-?,Km/{2\f{zo)\ ^^..^ v{zo) + 3Km/{2\f{zo)\ 

In particular, in view of ([5^ 

V^(^o) - 3/4 ^ ^ ^ vjzo) + 3/4 

Furthermore, L is a continuous function on the domain {z G C | \f{z)\ > 
max(e, 2mK)}. 



Remark 5.1 Here we deal mainly with the asymptotics of entire solutions of ^5^) 
However some statements are valid for arbitrary solutions of ^5^). In particular, for 
validity (E^j and no assumptions on the smoothness are needed. 
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5.6. Asymptotics and dynamics in the real case. Properties of the Julia set 

Further refinements are possible when A > 1 is real and P{z) = Pmz"^ + ■ ■ ■ + piz + po 
is a polynomial with real coefficients. 

Without loss of generality we can assume also that: 

/(O) = P(0) = 0; P'(0) = A > 1 and /'(O) = 1 

For reader convenience, we recall now some basic notions from the iteration theory 
and complex dynamics |25 |I110] . 

To make things shorter, we give here definitions which slightly differ from standard 
ones, but are equivalent to them in the polynomial case. That will suffice our needs, in 
what follows. 

We will especially need the component of oo of Fatou set J-'{P) given by 

J^oo(P) = G C I lim P(")(;z) = oo] , (61) 

The filled Julia set is given by 

IC{p) = {zeC \ (p(")(^))„eN is bounded} = C\T^{P). (62) 

Now, according to jllOj one can define the Julia set J7'(p), as 

dJCiP) = dT^iP) = J{P). (63) 

In the case of polynomials this can be used as an equivalent definition of the Julia set. 

Let f{z) be an entire solution of fl50|) . In contrast with the previous section (where 
asymptotics of 1/(^)1 is studied) here we collect some results on the asymptotics of the 
solution f{z) itself (in some angular regions of the complex plane). Our presentation is 
based on ( [MIIHIIES] ) . 

Theorem 5.4 ( |69L Theorem 1]) Let f he an entire solution of the functional 
equation ( fJOl) . Furthermore, suppose thatj^oo, the Fatou component of oo of P , contains 
an angular region of the form 

Wp = {zeC\{Q}\ larg^l </3} 

for some /3 > 0. Then for any e > and any M > the asymptotic relation 

/W=exp(.^g({^)+o(|.r^O) (64) 

holds uniformly for z G Wp^^! where Q is a periodic holomorphic function of period 1 
on the strip {w E C \ \^w\ < j^}. The real part of z^Q{^^) is bounded between two 
positive constants; Q takes real values on the real axis. 

Remark 5.2 Notice that the condition on the Fatou component J^^o is used in the proof 
of this theorem to ensure that f{z) tends to infinity in the angular region Wp. Therefore, 
this condition could be replaced by 

lim f{z) = oo for \ argz| < /3. 

z— >oo 
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Yet a stronger result can be derived under the additional assumption that the Julia 
set J{P) of polynomial P{z) is real (see Corollary 15. below). 

Remark 5.3 The latter assumption on the reality of J{P) looks artificial at the first 
glance, hut, in fact is typical in applications, related to the diffusion on fractals. The 
rough explanation of the latter fact is the following: 

The zeros of the solution f{z) for [5U) are eigenvalues of the infinitesimal generator 
of the diffusion, i.e. the "Laplacian on the fractal". The latter operator is self-adjoint 
, and its eigenvalues are real. Therefore zeros of f{z) are real and, finally, this implies 
that Julia set of polynomial P{z) should he real. This motivates our special interest for 
Poincare equations with polynomial P{z), having real Julia set. 

Corollary 5.1 ( \S7\ Corollary 4.1]) Assume that P is a real polynomial such that 
J{P) is real and all coefficients pi (i > 2) of P are non-negative. Then J{P) C M~U{0} 
and therefore 

f[z)^e^v(z''Qi^^\] (65) 



for z — 7- oo and \ argz| < vr. Here Q is a periodic function of period 1 holomorphic in 
the strip given hy < j^^. Furthermore, for every e > the real part o/ 2;''Q(|^|^) 
is hounded hetween two positive constants for \ argz| < vr — e. 

If P{z) is a quadratic polynomial ( a case, arising in numerous applications) it is 
possible to give an exact criterion for reality of J{P)'. 

Lemma 5.1 ( |95|, Lemma 6.7]) Let 

P{z) = az{z -u), Oy^ueR (66) 
Then Julia set J{P) is real, if and only if the following condition is fulfilled 

I I ^ / 2> ^ > 

a\uj\ > < 67 
' ' - [ 4, w < ^ ^ 

It would be interesting to find some constructive conditions for reality of J{P) in 
terms of coefficients of P (at least for polynomial of small order in the beginning). 
Note also, that in the above Lemma |-P'(0)| = a\u:\ and m = deg P{z) = 2. 
Therefore (1671) can be rewritten in the form: 



inO)l> (68) 

It turns out that the necessity part of latter criterion (in the form ( 168|) ) is valid 
for general polynomial P, with real Julia set J{P). Namely, the following result of 
Pommerenke-Levin-Eremenko-Yoccoz type (on inequalities for multipliers ) is true: 
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Theorem 5.5 ( |87L Theorem 4.1]) Let P be a polynomial of degree m > 1 with 
real Julia set J{P)- Then for any fixed point C, of P with min J7'(P) < ^ < max J'{P) 
we have \P'{^)\ > m. Furthermore, |P'(min J7'(P))| > and \P'{maxJ'{P))\ > m? . 
Equality in one of these inequalities implies that P is linearly conjugate to the Chebyshev 
polynomial Tm of degree m. 

Remark 5.4 This theorem can be compared to 11 05[\ 1 09^ 111^11^ where inequalities (of 
the opposite direction ) for the multipliers of P with connected Julia sets were derived. 

In the foregoing we deah mainly with the asymptotics of solutions for the polynomial 
Poincare equation, only. 

The study of the asymptotic behaviour of meromorphic solutions of the general 
Poincare equation 

fiXz) = Rifiz)), zee, 

where R{z) is rational function (rather than polynomial P{z)) still, to large extent, 
remains an open challenge. 
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